Oscillatory pairing of fermions in spin-split traps 
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As a means of realizing oscillatory pairing between fermions, we study superfiuid pairing be- 
tween two fermion "spin" species that are confined to adjustable spin-dependent trapping poten- 
tials. Focusing on the one-dimensional limit, we find that with increasing separation between the 
spin-dependent traps the fermions exhibit distinct phases, including a fully paired phase, a spin- 
imbalanced phase with oscillatory pairing, and an unpaired fully spin-polarized phase. We obtain 
the phase diagram of fermions in such a spin-split trap and discuss signatures of these phases in 
cold-atom experiments. 
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I. INTRODUCTION 

The idea that Cooper pairing in the presence of a den- 
sity imbalance of two interacting fermion species natu- 
rally yields oscillatory pairing correlations in real space 
was put forth decades ago. However, to date, this phe- 
nomenon, known as Fulde-Ferrcll-Larkin-Ovchinnikov 
(FFLO) pairing [l], has not been conclusively ob- 
served. (Related effects have been clearly seen in 
superconductor-ferromagnet hybrid systems where the 
proximity-induced pair correlations in the ferromagnet 
exhibit oscillations 0].) In recent years, atomic physics 
experiments have explored paired fermion superfluidity 
in cold atomic gases jJ-S] , a new setting for the observa- 
tion of FFLO pairing correlations under a density imbal- 
ance between the two "spin" species — a possibility that 
has inspired a large amount of recent theoretical and ex- 
perimental activity Q ■ Much of the excitement follows 
from the extreme tunability of cold-fcrmion experiments, 
which exhibit several experimentally-adjustable parame- 
ters including the interactions, the densities of the dif- 
ferent species, and the trap geometry. Of late, atten- 
tion has focused on one-dimensional (ID) systems with 
global spin imbalance or spin-dependent poten- 

tials (l5lll6j , where the parameter regime occupied by the 
FFLO state is predicted to be significantly wider than in 
the three-dimensional (3D) case Indeed, recent 

experiments (2fj| on quasi-lD spin-imbalanced fermionic 
gases have observed a partially polarized state, although 
associated oscillatory pairing correlations have yet to be 
confirmed. 

In this article, we propose a new ID setup to achieve 
FFLO pairing in cold atomic gases: a balanced mix- 
ture of two hyperfine species of attractively interacting 
fermionic atoms that are separately trapped in a control- 
lable way, as illustrated in Fig. [IJa) — a situation we call 
a "spin-split trap." This setup provides an effective spa- 
tially varying chemical potential difference between the 
two spin states due to the separate trapping potentials 
and yields an alternate, dynamically controllable route to 
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FIG. 1: (Color online) (a) Illustration of our proposed spin- 
split trap setup showing separate trapping potentials for two 
fermion species, t and L (b) Numerically-determined spatial 
profile of the pairing gap A(z) (solid blue curve, axis on left- 
hand side of graph), total density p(z), and spin imbalance 
M(z) (dashed light-brown and solid red curves, respectively, 
axis on right-hand side of graph, normalized by max z [p]), 
showing oscillatory pairing along with a local imbalance. 



achieving oscillatory FFLO-like pair correlations in cold 
atomic gases, controlled not by an imposed global popu- 
lation imbalance but, rather, by the separation between 
the two traps and the ensuing local imbalance. 

The spin-split trap, whose 3D counterpart was studied 
in Rcf . [2l| , is described by the spin-dependent potentials 
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V a {z) = -mui z (z — aoiy 



(1) 



where u) z is the trapping frequency, m is the atomic mass 
and a = ± corresponds to the two hyperfine species. 
Thus, the centers of the two traps are separated by a 
distance 2d. For d — > 0, the ground state is a singlet 
s-wave superfiuid with a vanishing spin imbalance every- 
where in the cloud. As argued below using local density 
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FIG. 2: (Color online) The local properties of the system in 
the spin-split trap can be understood using the phase dia- 
gram of the uniform imbalanced system, taken from Ref. [{J, 
showing fully paired (FP), partially polarized (PP), and fully 
polarized (FPo) phases as well as the vacuum. (Here p and 
h are measured in units of mg 2 /ih 2 , where g is the ID cou- 
pling constant.) The red curves A, B, and C represent the 
LDA trajectories followed as a function of z by the spin-split 
system for d < d c , d = d c , and d > d c , respectively. 
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FIG. 3: (Color online) (a)-(d) Spatial profile of the gap, total 
density, and magnetization in the z > region [represented as 
in Fig.GJb)] at d = 0, 0.176, 0.182, and 0.9R Z , respectively. In 
addition, in (a), the gap function obtained by BCS plus LDA 
is shown (dashed blue curve). In (b) and (c), the separations 
are just below and above the critical value d c for appearance 
of the first node. Note the scale change on the z axis in (d) . 



arguments and a Bogoliubov-de Gennes (BdG) treat- 
ment, for nonzero d, however, the split traps promote 
a local spin imbalance. We find that beyond a critical 
separation, d > d Cl the split-trap geometry displays os- 
cillatory pairing correlations, as depicted in Fig. [TJb), 
which shows the local pairing amplitude A(z), total den- 
sity p(z) = p^{z) + Pi(z), and magnetization (spin im- 
balance) M(z) = pf(z) — Pi(z). 

II. LOCAL DENSITY APPROXIMATION 

An intuitive understanding of the spin-split-trap sys- 
tem can be found using the local density approxima- 
tion (LDA) along with the known behavior of the ho- 
mogeneou s sp in-imbalanced gas derived using the Bethe 
ansatz 

HE!- The phase dia gram, shown in Fig. [5J dis- 
plays three phases as a function of the net chemical po- 
tential p = (/i-j- + Pi)/2 versus the chemical potential im- 
balance (magnetic field) h = (/x-j- — p^)/2, namely a fully 
paired (FP) state, a fully polarized (FPo) state, and a 
partially polarized (PP) state. The PP state is expected 
to be of the FFLO type, having an oscillatory pairing 
amplitude [1, [TU , as corroborated by our studies below. 

Within LDA, the trapping potential in our system en- 
ters as a spin-dependent spatially-varying chemical po- 
tential, Pa-(z) = po — V a (z), where po is the global 
chemical potential of the system. For the harmonic 
trap of Eq. ([1]), p and h are then related through 
p = po — h 2 1 (2mLu z d 2 ) , which corresponds to downward- 
facing parabolas in the p versus h phase diagram. In 
Fig. [2] we show three curves corresponding to different 
values of the separation d. One can see that they traverse 
different phases from the center z = (where h = 0) to 
the edges of the trap. For small separation d, the sys- 



tem is described by a tight parabola and is thus con- 
fined to the fully paired phase, but with increasing d, the 
parabola broadens and, beyond a critical separation d Cl 
traverses all three phases as a function of position. In this 
case, at small z, the local potential imbalance h remains 
small enough that the system is (locally) fully paired. At 
larger z, the local h exceeds a critical value such that 
(locally) the system enters the PP phase. At even larger 
z, near the edges of the trap, the system is (locally) in 
a fully polarized normal phase. Thus, the system con- 
currently hosts all three phases. Note that, in contrast, 
in the case of a globally spin-imbalanced system with a 
single trap, the system traces a vertical line in the phase 
diagram, yielding two regions — a partially polarized core 
and either fully polarized or fully paired edges @. 



III. MICROSCOPIC THEORY 

We now model the spin-split system using a micro- 
scopic description which enables a more detailed anal- 
ysis, confirms the salient features described above, and 
shows a direct correspondence between local spin imbal- 
ance and oscillatory pairing. We study two species of 
fermions, ip^^(z), in a ID harmonic potential charac- 
terized by the trapping frequency lo z . In atomic sys- 
tems, this limit can be achieved in a highly anisotropic 
trap with a transverse trapping fre que ncy uj r such that 
Nu z /uj r < 1 and N\a s \/R z < 1 0,122]. Here, AT is t he 
number of fermions of each spin species, R z = \f2N— 1 £ z 
(with £ z = yjh/muj z the oscillator length) is the classical 
radius of the free gas in the z direction, and a s is the 
s-wave scattering length for the two-body interactions. 
The system is then described by the effective ID Hamil- 
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FIG. 4: (Color online)Phase diagram as a function of separa- 
tion and temperature (N = 40, g/huj z R z — 1). The solid line 
separates the normal phase and the superfluid phase. In the 
superfluid phase, the dashed lines separate regions of the gap 
functions with different number of nodes. 

tonian, 



(2) 



where H® = — (h 2 /2m)d 2 + V a (z) — [lo is the one-particle 
Hamiltonian. The ID coupling constant is given as g — 
2h 2 a s /[m£ 2 (l — l.Q33a s /£ r )] with the transverse oscillator 
length l r = y/h/mcor [H]. 

We analyze our system within the standard BdG treat- 
ment, which has been widely applied to the imbalanced 
system [24| , taking into account spin-dependent trapping. 
The mean-field Hamiltonian, which self-consistently in- 
corporates the Hartree potential U a = g^>\^) a ) and pair- 
ing gap A = g(ipiipt), takes the form 

H M = [ dz[Y,^HH° a + U a )4>„ + (A^l + K.c.)]. (3) 

We obtain the extended BdG equations in the quasi par- 
ticle eigenbasis by a spin-dependent Bogoliubov trans- 
formation, tpa(z) = J2l u ncr(z)%a ~ ^no-C^Tn.-o-]- Wc 
n 

use an iterative numerical procedure [25[ to find self- 
consistent solutions for p a {z) and A(z). Parity symme- 
try between the potentials of the two species, Vj. (z) = 
V\ (— z), ensures parity symmetry of the gap function; we 
find that the even-parity solution, A(z) = A(— z), is al- 
ways energetically favorable. The data presented in the 
following were obtained for N — 40 and g/ftw z R z = 1. 



IV. RESULTS 

We first focus on the manner in which oscillatory pair- 
ing correlations emerge with increasing separation d. In 
Fig. [3l we show the pairing gap A(z), total density p{z), 
and magnetization M(z) for a sequence of four spin- 
split-trap systems with increasing d. Panel (a) shows 
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FIG. 5: (Color online) (a) Density plot of the polarization P 
as a function of position z and separation d. The gray scale 
is bounded by and 1. The dashed white (black) contours 
correspond to P = 0.01 (0.99). The solid curves indicate 
positions of the first four nodes, (b) Momentum distribution 
n(k) for d = (solid blue curve) and 0.25i? z (dashed red 
curve), normalized by n(0) at d = 0. (c) Fraction of pairs fo 
within the central peak [see (b)] of n(k) (solid curve, axis on 
left-hand side of graph) and positions of the first four nodes 
(dashed curves, axis on right-hand side of graph) vs d. 



the d = case which is fully paired with M = ev- 
erywhere, as expected. The non-monotonicity of A(z) 
roughly reflects the functional dependence of the ID 
BCS gap on the l ocal chemical pot ential /i, that is, 
A(z) oc fi(z) cxp[— [i{z) I (mg 2 )]. Panel (b) shows 

that a small separation, d < d c , does not lead to qual- 
itative changes of the pairing correlations and the mag- 
netization. Here, the local h is small enough everywhere 
that it is energetically favorable for the system to remain 
fully paired (i.e., the system is below the Clogston limit). 
Panel (c) shows the system just beyond the critical sepa- 
ration d Cl such that, near the edge of the cloud where the 
local h is largest and of order A, the gap function A(z) 
exhibits a node and the magnetization is finite. As d 
increases further, the region of oscillatory FFLO correla- 
tions increases and more nodes appear. The progression 
of nodes is captured in Figs. [TJb), O^c), and Eld). Ini- 
tially the number of nodes increases as d increases, but 
then, beyond a characteristic distance of the order of the 
cloud size, diminishes before the system fully separates 
and becomes normal. 

We find that the nodal structure is robust against fi- 
nite temperature effects. This is illustrated in the global 
phase diagram in Fig. |4l obtained using the parameter 
values specified above. Within the superfluid phase, re- 
gions with different numbers of nodes in A(z) are in- 
dicated. We note that the transition temperature in 
the spatially modulated phase is of the same order as 
in the fully paired phase. The number of nodes de- 
creases with increase in temperature, consistent with the 
shrinking of the FFLO region in globally imbalanced sys- 
tems [26|. As for trends with variation of the system 
parameters, we numerically find that the critical separa- 
tion d c = d c j R z is independent of N and linearly depen- 
dent on g = g/(hu) z R z ) around 5 = 1 (in the regime of 
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numerical convergence), which is consistent with rough 
estimates based on BCS combined with LDA. 

Our results for the behavior of interacting fermions 
in the spin-split trap clearly show the intimate connec- 
tion between a nonzero polarization and oscillatory pair- 
ing correlations. In Fig. [5ja) we show the polarization, 
P(z) = M(z)/p(z), as a function of position, z, and sep- 
aration, d, along with the spatial position of the nodes 
in A(z). It can be seen that the nodes exist only in the 
partially polarized region, < P < 1. The correlation 
between the polarization and nodal structure indicates 
that this region is indeed of the FFLO type, and is sur- 
rounded by a fully gapped superfluid for P — > toward 
the center of the spin-split trap and a fully polarized nor- 
mal fluid for P — > 1 at the edges. 

V. EXPERIMENTAL ASPECTS 

A direct measure of oscillatory pairing is the pair mo- 
mentum distribution function defined as 

n(k) = J dzdz'e^ z - z '\^(z)^l(z)^(z')^(z')), (4) 

which is experimentally measurable in dynamic- 
projection experiments [231. In the homogeneous case, 
the FFLO phase is characterized by a peak in n(k) at a 
characteristic nonzero wave vector k that depends on the 
spin imbalance [TlLIT3|. Typical plots of n(k) in the spin- 
split trap are shown in Fig. [5jb) for the cases of d = 
and d > d c . Due to the spatial inhomogencity of the po- 
tential imbalance h, the system does not possess a char- 
acteristic wavevector. However, n{k) undergoes sudden 
changes with increasing separation as Cooper pairs are 
shifted to higher momenta. As shown in Fig. El[c), the 
weight under the central peak suddenly decreases each 
time a new node appears in A(z). Thus, n(k) displays a 
striking signature of the modulated phase. 

We now turn to the issue of experimentally realizing 
a spin-split-trap system. This setup can be achieved via 
spin- selective trapping potentials [5^, [29| . Additionally, 
a tunable spin-split trap may be achieved using a mag- 
netic field gradient [22[ l3(il l3l| , exploiting the distinct 
hypcrfine-Zeeman states of the two fermion species. To 
see this, we note that the competition between the Zee- 
man effect and hypcrfine interaction leads to a nonlin- 
ear energy difference between the two spin states mp±. 
We use the Breit-Rabi formula [32| to find the spatially- 
varying part of the energy difference AV(z) = V\(z) — 
Vl(z) in the presence of a field gradient. Assuming a 
spatially-varying field of the form B(z) = B + B'z and 
expanding the Breit-Rabi formula near the background 
field B, we obtain AV(z) = 2ttHB' fl(B)z with p, the ef- 
fective "magnetic moment" given by 

2mFg I _B_ 

m = f e * i T ;° n2 ■ (5) 

2 — ' /i , 4mjy B , B 1 
<y=± y ' 2/+1 B + "Bj 



Here, /xb is the Bohr magneton, / is the nuclear spin, Bq 
is the hypcrfine field, and g ~ 2. 

Using Eq. ([T]), we see that a spatial separation d re- 
quires a field gradient B' = 2mui^d/p,. In the case of 
interest, we expect that B is close to a Fcshbach reso- 
nance (FR) in order to enhance T c and that B' is small 
enough that o s can be treated spatially independent in 
the system. For 6 Li (/ = 1, Bp = 81 G), using the hyper- 
fine levels mp± = —§(+5) [33| near the FR at B = 691 
G, we find /in ~ 6 x lO~ 3 [j, B . Assuming a typical trap 
frequency u> z ~ 2ir x 100 Hz, a field gradient B' Li of the 
order of a few hundred G/cm can achieve a separation 
d of a few £ z . (The required field gradient for the more 
commonly used two lowest hyperfine levels of 6 Li is about 
an order of magnitude larger and thus much less experi- 
mentally viable.) The most promising case is that of 40 K 
(I = 4, B — -459 G) with m F ± = ~|(— §) near the 
FR at B = 202 G. Here /2k — 0.1/Ub and, for the same 
uj z as above, the required gradient B' K ~ 50 G/cm. 



VI. CONCLUSION 

In summary, we have proposed a novel setting, the 
spin-split trap, for observing FFLO-like oscillatory pair- 
ing correlations, driven by a local density imbalance due 
to the separate trapping potentials of the two fermion 
species. Our BdG calculations, supported by LDA, show 
that the competition between the tendency to pair and 
the tendency towards forming a spin imbalance leads to 
a rich structure that is revealed in quantities such as 
the local pairing amplitude and magnetization, as well 
as in the pair momentum distribution. Immediate fu- 
ture directions include investigating the spin-split system 
through other techniques amenable to ID, such as den- 
sity matrix renormalization group (DMRG) and quan- 
tum Monte Carlo methods, and exploring the exciting 
prospect of coupling arrays of spin-split ID systems. 
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